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Abstract 

The dynamics of a central spin- 1/2 in presence of a local magnetic field and a bath of 
N spin- 1/2 particles is studied in the thermodynamic limit. The interaction between the 
spins is Heisenberg XY type and the bath is considered to be a perfect thermal reservoir. 
In this case, the evolution of the populations of the reduced density matrix are obtained 
for different temperatures. A Born approximation is made but not a Markov approximation 
resulting a non-Markovian dynamics. The measure of the way that the system mixes is 
obtained by means of the von Neumann entropy. For low temperatures, results show that 
there are oscillations of populations and of the von Neumann entropy, indicating that the 
central spin becomes a pure state with characteristic time periods in which it is possible to 
extract or recuperate information. In the regime of high temperatures, the evolution shows 
a final maximum mixed state with entropy S = In 2 as it is expected for a two level system. 

I Introduction 

The most promising systems that could be scaled to make practical realizations of quantum 
computation and quantum information are the spin systems in solid state nanostructures [lj 
[21 [3l S]. Developments in the last few recent years, permit to manipulate, tune and have full 
control over the spin of individual defined number of electrons in semiconductor nanostructures 
like quantum dota3 [5j [6]. This particular experimental abilities of control over the spin had 
motivated the study of individual spins as the fundamental system [TJ. However, the correlations 
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of the spin with many degrees of freedom of the surrounding environment lead to finite lifetime 
of quantum superpositions, bringing pure states into mixed ones. 

The principal source of noise in this solid state spin nanodevices arises from the hyperfine 
interaction with nuclear spins; a lot of efforts has been devoted to model spin bath systems 
El HO]. Different theoretical proposals have been made in order to model these systems, 
however, a special interest have been devoted to the configuration known as spin star network 
[111 I12j . This configuration consists of a N + 1 spin- 1/2 particles, where the system of interest 
is a central spin and the remaining N spins surround the central spin at equal distances. This 
configuration implies a spatial symmetry which allows a formal and exact analytical solution. 
This N surrounding spins act as a thermal bath. 

In the same way, different theoretical calculations have been made in order to obtain the reduced 
dynamics of the central spin. Many of these works have considered Ising type interaction [13|, I14j . 
The work of reference [TT], consider a star network where the interaction between the spins of 
the bath is neglected. A similar situation is considered in reference [12], where the interaction 
between the central spin and the bath is Heisenberg XY type and the the bath is taken in an 
unpolarized infinite temperature state which in practice is not a feasible realization because 
of the Coulomb blockade [7J. In the reference [15] it is considered a more general case of the 
dynamics of two central interacting spins in a bath of interacting spins in order to describe 
the entanglement between the two central spins without considering the entanglement with the 
bath. 

In this paper, we consider the reduced dynamics of one central spin in a bath of interacting spins 
in a configuration of spin star network. The central spin interacts via Zeeman effect with a local 
magnetic field, the interaction between spins is Heisenberg XY type and the thermodynamic limit 
is considered. The methodology is similar to that of [15], however, we consider the dynamics 
of one single spin and the way that it becomes mixed. In addition, in reference [15] it is 
considered a spin wave theory by means of the Holstein-Primakoff transformation |17j and we 
study the problem of small fluctuations in the context of Schwinger bosons [15] . The Schwinger 
transformation and the Holstein-Primakoff are closely related as it is shown, however, Schwinger 
bosons have more natural physical interpretation. The principal reason to consider the system 
of one central spin is that they are really controllable systems [HI El E] and fulfill the criteria of 
DiVincenzo for universal quantum computing [U [19] . 

The paper is organized as follows. In Section[H]we present the model by means of the Schwinger's 
oscillator model of angular momentum [18] and in the thermodynamic limit we obtain a Jaynes- 
Cummings Hamiltonian type. Then the reduced dynamics of the central spin and the evolution 
of the populations of the central spin is obtained. In Section [nil we calculate the time evolution 
of the von Newmann entropy. Different temperature regimes are considered showing, how the 
temperature influences the mixing of the state. Conclusions are drawn in Section [TVl 
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II The model 



The characteristics of the model are similar to that of [T5], but our principal interest is to 
describe the coupling of the central spin with the environment. We consider a spin star network 
of N + 1 interacting spin-1/2 particles. The interaction is Heisenberg XY type and the system 
of interest is the central spin in presence of a local magnetic field. The bath is represented by 
the N 1/2 spin particles surrounding the central spin. The total Hamiltonian is 

H = H S + H B + H SB , (1) 

where Hs and Hb are the Hamiltonians of the system and the bath respectively. The interaction 
Hamiltonian between the system and the bath is denoted by H$b- Reminding that there is a 
local magnetic field interacting only with the central spin, and that all particles interact via 
Heisenberg XY, each of this terms can be written as [15] 



H S = fioS', (2) 

^ = #E^ +5 7 + s r^;> (3) 

and 

r n JV 
Hsb = -% S+J>- + 5 ~E^ + , (4) 

v iV L i=l i=l 

where, /j,q represents the local magnetic field in the z direction, g is the coupling constant between 
the spins of the bath, and go is the coupling constant between the central spin and the bath spin 
particles. We should stress that the coupling constants has been scaled as g/N and g^/^/N in 
order to obtain the thermodynamic limit [151 120]. The correspondig raising and lowering spin 
operators, written in terms of the usual spin operators S x , S y , are = S x ± iS y . The total 
spin of the bath can be written as J± = YliLi ^ > so that the central spin couples to an effective 
collective bath of angular momentum J. In terms of J± the bath and the interaction terms of 
the Hamiltonian become: 

H B = jj(J + J- + J-J + )-g (5) 

and 

H SB = -^{S + J^ + S-J + ) (6) 



respectively. The spin-i operators can be mapped to Bose type operators by means of the 
connection between the algebra of the angular momentum and the algebra of two independent 
oscillators. The problem has been studied by J. Schwinger and the process corresponds to spin 
wave theory in order to describe the spins in terms of small fluctuations [18j. 
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II. 1 Schwinger's bosons 

Consider two non-coupled harmonic oscillators, or Schwinger bosons, denoted by b and a. The 
usual number operators are given by 

N b = tfb, N a = a^a. (7) 

It is considered that b and a are independent and follow the Bose statistics: 

[6,6t] = l [a,at] = l, (8) 

[b, a f ] = [a, 6 f ] = 0. (9) 

Then iVj, and N a conmute, [N b , N a ] = 0, and there exist simultaneous eigenkets | n b ,n a ) where 
the correspondig actions of the b and a operators are as ussually: 

b 1 | n b ,n a ) = sfn h + 1 | n b + l,n a ), a f | n 6 , n a ) = V™a + 1 I ™b, n a + 1) (10) 

and 

b | n 6 , n a ) = jn~b \n b -\, n a ), a \ n b , n a ) = jn~ a \ n b , n a - 1). (11) 

By successive applications of the creation operators and a) we can obtain the most general 
eigenkets, 



(&t)»( 



\n b ,n a )= " ' K ' | 0,0). (12) 
The Schwinger transformation is defined as (with ft = 1) 

J+ = 6 f a J_ = a f 6. (13) 

and 

J z = \(N b -N a ). (14) 

It is easy to verify that while b and a follows a Bose statistics, the J± satisfy the usual relations 
of commutation of angular momentum, 

[J+,J_] = 2J Z , (15) 

[J Z1 J±]=±J± (16) 
and acting on the common base | n b ,n a ) they give 



J+ | n b , n a ) = y/(n b + l)n a \ n b + 1, n a - 1), (17) 
J_ | n&, n ) = V n 6 (n a + 1) | n b - 1, re a ), (18) 
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J z I n b ,n a ) = -(n b -n a ) \ n b ,n a ). 



(19) 



The Schwinger transformation has an immediate physical interpretation. There are two kinds 
of collective excitations, and equation (|19p suggests that n b is the number of spins up, while n a 
is the number of spins down. So, the two kind of collective excitations are: creating units of 
spin up and creating units of spin down. The meaning of J + is that it destroys one unit of 
spin down at the same time that it creates a unit of spin up. Similarly, J_ destroys a unit of spin 
up and creates one of spin down. The different processes occur without interactions between 
the two excitations as has been established by equation Q. 
On the other hand, the spin magnitude defines the physical subspace 




(20) 



and the Schwinger transformation can be written as [16] 



J+ = bWN-tfb, 




J_ = - tfb b. 




The Hamiltonian of the bath in terms of the Schwinger bosons can be written as 




(23) 



Since N = b^b + a) a and a) a = N — N b , we can write the bath Hamiltonian as 




(24) 



In the thermodynamic limit, N — > oo, the Hamiltonian of the bath becomes 



H B = 2gtfb. 



(25) 



The interaction Hamiltonian H$b given by equation (J6j) becomes 




(26) 



and using equations (|21|) and (|22|) we obtain 




(27) 



which in the thermodynamic limit, N — > oo gives 




(28) 
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and the total Hamiltonian can be written as 



H = p S z + 2gtfb + g (S + b + S~ 6 f ). (29) 

This Hamiltonian is equivalent to the Hamiltonian of a two level atom in presence of a bath 
of bosons of one single mode. It is clear that the Hamiltonian (J29]) is analog to the quantum 
model of radiation-matter interaction with a single mode in the rotating wave approximation: 
the Jaynes-Cummings model. 

The problem of spin wave theory can be studied by means of the Holstein-Primakoff transfor- 
mation instead of Schwinger bosons. In this case the angular momentum operators are mapped 
to boson operators by [TT] 

J+ = tfy/2S-N b , (30) 
J_ = V25 - N b b, (31) 

and 

J z = N - N a . (32) 
It is easy to show that the operators J obeys the angular momentum commutation relations, 

[J a ,J f3 ]=ie af3 ''J^ (33) 

The two transformations known as Schwinger bosons, SB, and Holstein Primakoff, HP, are 
related by the following correspondence 

SB -> HP 
b b 
a -» VN - N b . 

As was discused before, this is a consequence of the physical subspace spanned. However, 
the Schwinger bosons transformation has a natural physical interpretation, as it was presented 
above. 



II. 2 Reduced dynamics 

Let us consider the solution of the dynamics of the density operator of the system to obtain 
the dynamics of the central spin. Suppose that the total system represented by pit) at t = is 
separable, p(0) = \tp)(i^\ <8> ps, where the initial state of the central spin is in a pure state, \tp), 
and the bath is in a thermal bath given by the density operator pxh- 

exp(-H B )/T 

PB = PTh = ^ , (34) 
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where the Boltzmann constant has been taken equal to one. The Born approximation is made in 
the sense that the state of the bath is time independent, pB(t) = Pb for any time. The partition 
function Z is given by 

Z = Tr[exp(-H B /T)}, (35) 

where the Hamiltonian of the bath, after the Schwinger transformation, is given by Hb = 2gb^b. 
Then the partition function can be calculated analytically to obtain 

The dynamics of the total density operator of the system can be obtained from the Liouville-von 
Neumann equation, 

^ = -i[H,p), (37) 

where H is the total Hamiltonian given by the equation (|29p . This Hamiltonian can be written 

as 

H = H 1 + H 2 (38) 

where 

H 1 = 2gtfb + 2gS z (39) 

H 2 = AS* + g (S + b + S- &t). (40) 

where A = p Q — 2g, usually known in quantum optics as detuning. The two parts of this 
Hamiltonian commute, [Hi, H 2 ] = 0, and the time evolution operator can be written as 

U(t) = exp(-iHt) = exp(-iffii) exp(-iH 2 t) = U^U^t). (41) 

In the two dimensional subspace spanned by the eigenvectors of S z , that is in the base {| |), | J,)}, 
where | f) and | j) indicate spin up and spin down respectively, the first factor, U\(t), of the 
time evolution operator is diagonal. On the other hand, the second factor, [/^(t), is expanded 
in order to obtain, 

/ cos(^vWTI) -ib sH9 J^f b) \ 

within the limit of zero detuning, po = 2g, that is, with a fixed magnetic field of magnitude twice 
the magnitude of the internal interaction. In the two dimensional base,{||), the evolution 
of the total density operator is obtained by 

p(t) = U 2 (t)p(p)U}(t), (43) 
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and the central spin dynamics is obtained by tracing over the degrees of freedom of the environ- 
ment, ps = Tvb{p) [21] ■ The upper level population is given by 

oo 

Puu{t) = J2 Pn cos2 (9o Wn + 1), (44) 

n=0 

while pm = 1 - p uu , with p uu = (| \ps(t)\ |) and p dd = {[ \ps(t)\ j). The coefficients P n , are 
the probability of the thermal mode P n = (n\prh\n)- In this case 

Note that the coefficients of the populations P n include the effects of the interaction between 
the bath spins and the thermal energy T. Figure Q] shows the temporal evolution of the density 




Figure 1: The upper left plot shows the time evolution of the element p uu for three different 
temperatures T = 40g, T = lOg and T = O.lg with time in units of g^ 1 . The initial state of 
the central spin is ||). The upper right plot shows the same result for a larger time scale and 
for the temperatures of T = 40g and T = Wg. The down left plot shows the result obtained for 
a temperature of T = lOOOg, with a faster relaxation process. The lower right plot shows the 
graph but fora larger period of time and for a temperature of T = lOg. 

matrix element, p U u = {\ \ps{t)\ T)- The time is in units of g^ 1 and the temperature in units 
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of g. The initial state of the system is spin up ||). The effect of the temperature is clear 
in the sense that for increasing temperature the relaxation process is faster. The plots show 
clear oscillations of the relaxation process which get smaller as temperature increases. This is 
recognized as evidence of non Markovian dynamics (for example see reference [22]). The period 
of oscillation is approximately 3 units of g^ 1 for almost all the temperatures and its amplitude 
decreases as the temperature increases, showing that the general effect of the temperature is to 
bring the system from pure states to mixed states. Different temperatures have been considered. 
The upper left plot shows the time evolution of the element p uu for three different temperatures, 
T = 40<7, T = lOg and T = O.lg with time in units of g^ 1 . The graph of the right shows the 
revivals of dynamics as a signal of the interchange of information between the system and the 
bath. The lower left plot shows the evolution of the population for a temperature of T = lOOOg. 
Although the temperature is high, in terms of the coupling constant g, it is possible to see that 
the memory effects still persist, but no so notable as in the case of low temperatures. The 
population oscillates around the value 0.5, for which the system becomes a completely mixed 
state. The lower left graph shows that the oscillations never disappear as it is clear from the 
result for the upper level population, equation (|44l) . 



Ill Von Neumann entropy dynamics 

The von Neumann entropy describes the departure of a system from a pure state or equivalently 
it measures the degree of mixture of a system. Taking two extreme values of zero for pure states 
and In TV for a maximally mixed state, N being the dimension of the Hilbert space. The von 
Neumann entropy is defined as 

S(p) = -Tr[plnp], (45) 

where p is the density matrix. An important property of the von Neumann entropy is that it is 
invariant under changes in the basis of p, and as we can see this property is primordial in order 
to calculate the evolution of the dunamics of the mixing. Since the reduced density matrix is a 
2x2 matrix, it can be written in terms of the Pauli matrices: 

where / is the identity matrix of dimension two, a are the traceless Pauli matrices, and s is 
known as the Bloch vector. The Bloch vector has magnitude one and its tip lies on the surface 
of the Bloch sphere when the state is pure, and it is easy to show that | s |< 1 in the case of 
mixed state. The eigenvalues Ai and A2 of the density matrix in terms of the components of the 
Bloch vector are: 

Ai = ^[l+|s|] 

2 (47) 

A 2 = i[i-|i|]. 
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On the other hand, from the reduced density matrix (|43p . we can write it as 

/ J2n=o p n cos 2 (gty/n + 1) % Y,n=o P n sin 2 (gty/n + 1) cos(gty/ n + 1) \ 

PS = 

\ ~i E^°=o p n sm 2 {gty/n + 1) cos(gty/n + 1) E^=o p n sin 2 (5ftVn + 1) / 

(48) 

where P n = e 2 9 n l T (\ — e 2f] l T ). Comparing the matrix given by equation (|46p with the matrix 
given by equation ([48]) we obtain the elements of the Bloch vector, the eigenvalues of p given by 
the equation (|47p . and the corresponding entropy of the system: 

S(t) = -A 1 ln(Ai)-A 2 ln(A 2 ). (49) 

In figure [2] we present the plot of the von Neumann entropy corresponding to the reduced 




Figure 2: Time evolution of the von Neumann entropy for the initial pure state ||) and for 
different values of the temperature, T = lg, T = 5g and T = lOg. The larger the temperature, 
the faster that entropy becomes S = In 2: the entropy of a two level mixed state. The system 
and the bath becomes nearly pure for low temperatures. 

density matrix for the central spin in the case of a initial pure state ||). As in figure [H the time 
is in units of g^ 1 , and the temperature in units of g. The result is notable in the sense that 
for low temperatures the central spin periodically returns to its initial pure state. In the case 
of higher temperatures the system goes more rapidly to a completely mixed state of entropy 
S = In 2 ~ 0.69, however, there remain oscillations with smaller amplitudes. In general, we can 
say that memory effects, or equivalently, non Markovian dynamics, are responsible of information 
gain for some particular time periods. The result obtained for the entropy in the regime of 
high temperatures is consistent with |12| . where the initial state was an infinite temperature 
unpolarized state. Some differences in the entropy of [12] occur due to different initial states for 
the central spin. 
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IV Conclusion 



The dynamics of a central spin in a bath of spins have been studied considering Heisenberg 
XY type interaction. An analytical solution was obtained for the reduced density matrix of 
the central spin in the thermodynamic limit. The bath acts as bosonic collective excitations of 
one single mode via Schwinger bosons. A separable state between the bath and the system was 
considered as initial state with the central spin in a pure state and the bath in a thermal state, 
evaluating the relaxation process for different temperatures. As the temperature increases, the 
relaxation process is faster as it is expected intuitively, however, memory or non Markovian 
effects are observed as oscillations of the populations which are reduced in amplitude with the 
increase of the temperature. The evolution of the degree of mixure of the system in presence of 
a bath was studied by means of the von Neumann entropy for different temperatures. At low 
temperatures the system becomes periodically completely pure. The increase of the temperature 
implies only partial mixing and leads the system to a complete mixed state with entropy S 1 = In 2 
in the case of infinity temperature. The results obtained coincide with that of the entropy of 
reference [12], where an unpolarized infinite temperature state was considered, but with the 
difference that in this case we can make calculations for different temperatures. It is important 
to note that a thermodynamic limit was considered here, however, as was shown in [12], the 
results do not depend on the number of particles, N, for N > 200 for longer time scales; and 
for smaller times scales for N ~ 20. We point out the advantage of the thermodynamic limit 
permitting analytical solution. 

We expect different behaviors of the entropy depending on the magnetic field as a control 
parameter of the degrre of mixing, as will be presented in a future work. On the other hand, 
our approach in terms of Schwinger bosons, permits to consider a general problem with clear 
physical interpretation giving also the possibility to describe a still more interesting system of 
interacting bath modes (work in progress). 
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